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Abstract. We present a complete library of synthetic spectra based on Kurucz's codes that covers the 2500-10 500 A wave- 
length range at resolving powers R P =20 000, 1 1 500 (= GAIA), 8500 (= RAVE), 2000 (= SLOAN) and uniform dispersions of 1 
and 10 A/pix. The library maps the whole HR diagram, exploring 51 288 combinations of atmospheric parameters spanning the 
ranges: 3500 < T cff < 47 500 K, 0.0< log£ < 5.0, -2.5 < [M/H] < 0.5, [a/Fe] = 0.0,-1-0.4, f =1,2,4 km sec" 1 , < V m < 500 km 
sec -1 . The spectra are available both as absolute fluxes as well as continuum normalized. Foreseen applications of the library 
are briefly discussed, including automatic classification of data from spectroscopic surveys (like RAVE, SLOAN, GAIA) and 
calibration of differential photometric indexes. Data distribution and access to the library via dedicated web page are outlined. 
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1. Introduction 

The interest in large and complete synthetic spectral libraries 
of normal stars ranges over many different application areas, 
including (a) automatic analysis and classification of large vol- 
umes of data, like those collected by ongoing spectral sur- 
veys, (b) derivation of radial velocities via cross-correlation 
against best matching templates, (c) calibration of spectro- 
scopic line/band classification criteria, and (d) calibration of 
photometric indexes. Tasks a-c require the synthetic spectra 
to match the absorption lines in observed spectra normalized 
to the local continuum, generally irrespective of how well the 
overall spectral energy distribution is reproduced. The goal of 
task d is instead to match the energy contained within pho- 
tometric bands, irrespective of the performance on the indi- 
vidual absorption lines. Therefore, synthetic spectral libraries 
well suited for one task does not necessarily perform equally 
well on another. For example, within the Kurucz's synthetic 
spectra suite, the inclusion in the computations of the so called 
'predicted lines' (hereafter PLs, Kurucz 1994) limits the range 
of applications. The PLs are calculated absorption lines due to 
transitions involving one and sometimes two levels whose loca- 
tions are predicted by atomic structure codes. Wavelengths for 
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The entire library of synthetic spectra is accessible via 
http://archives.pd.astro.it/2500-10500/ or (for the 1 A/pix version 
only) http://gaia.esa.int/spectralib/spectraliblA/ SpectraLibla.cfm 



these lines may be uncertain by >10 A, with large uncertain- 
ties affecting also their computed intensity. As a consequence 
the predicted lines may not correspond in position and size to 
the observable counterparts (cf. Bell et al. 1994, Castelli and 
Kurucz 2004). For all these reasons the predicted lines are not 
intended to be used for the analysis of high resolution spectra, 
but are an essential contribution to the total line blanketing in 
model atmospheres and spectral energy distribution computa- 
tions. Very accurate line wavelengths and shapes are in fact not 
needed for the statistical computations of the line opacity in 
model atmospheres as well as for the prediction of low resolu- 
tion spectra. 

Various libraries of synthetic spectra are becoming avail- 
able, and they are summarized and compared in their basic 
characteristics in Table 1 . These libraries widely differ in wave- 
length interval, spectral resolution, input model atmospheres 
and range of atmospheric parameters explored, so that their 
overlapping is not significant. 

As other libraries in Table 1, our one is based on Kurucz's 
atmospheres, line-lists and computing software. The main goal 
of our spectral library are spectroscopic applications of the a 
and b types above, which imposes for ex. the exclusion of the 
PLs from computations and therefore no pretension of accurate 
broad-band photometric performances particularly at the lower 
temperatures. Other relevant characteristics of our library are 
the various resolutions and rotational velocities at which it is 
provided, its wide wavelength range, the adoption of improved 
model atmospheres based on the new opacity distribution func- 
tions (ODFs) by Castelli and Kurucz (2003), the use of the TiO 
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Fig. 1. Coverage of the log g, Tgg plane by the synthetic spectral 
grid presented in this paper (figure available electronic only). 
Panel a: spectra with [a/Fe]=0.0 and £=2 km sec -1 ; Panel b: 
spectra with [a/Fe]=+0.4 and £=2 km sec -1 ; Panel c: spectra 
with [a/Fe]=+0.4 and f£2 km sec -1 

line list of Schwenke (1998), the inclusion of <?-element en- 
hancement and different micro-turbulent velocities. 

2. The spectra 

The whole grid of spectra in our library was com- 
puted using the SYNTHE code by Kurucz (Kurucz and 
Avrett 1981, Kurucz 1993), running under VMS operat- 
ing system on a Digital Alpha workstation in Asiago. 
We adopted as input model atmospheres the ODFNEW 
models (http://wwwuser.oat.ts.astro.it/castelli/grids/, Castelli 
and Kurucz 2003). They differ from the NOVER models 
(http://kurucz.harvard.edu/grids.html, Castelli et al. 1997) for 
the adoption of new ODFs, replacement of the solar abun- 
dances from Anders and Grevesse (1989) with those from 
Grevesse and Sauval (1998), and improvements in the molecu- 
lar opacities among which the adoption of the molecular line- 
lists of TiO by Schwenke (1998, as distributed by Kurucz 
1999a) and of H 2 by Partridge and Schwenke (1997, as 
distributed by Kurucz 1999b). For the combination of atmo- 
spheric parameters for which no ODFNEW models were avail- 
able at the time of writing, we adopted the corresponding 
NOVER input model atmospheres. We plan to update the li- 
brary from NOVER to ODFNEW models as more of the lat- 
ter will become available and to post the result on the li- 
brary web page (see below). We remark that both the NOVER 
and the ODFNEW models were computed with the over- 
shooting option for the mixing-length convection switched off, 
while Kurucz (K) atmospheric models (also available from 
http://kurucz.harvard.edu/grids.htmli were computed with the 
overshooting option switched on. Several papers have demon- 
strated that for stars with active convection (r e (j<9000 K) the 
no-overshooting convection treatment provides better agree- 
ment with the observations than the overshooting case does 
(Castelli et al. 1997, Smalley and Kupka 1997, Gardiner et al. 
1999, Smalley et al. 2002). The no-overshooting models used 
by us were computed for the mixing-length parameter to the 
scale height of 1.25. This value allows to fit the observed solar 
irradiance, while a lower value, like that of 0.5 suggested by 
Smalley et al., 2002, does it not. On the contrary, 0.5 seems to 
better fit the wings of H/3, provided that the position of the con- 
tinuum is known with an uncertainty smaller than 1%. In fact, 
a difference of 1 % in the position of the solar continuum cor- 
responds to the difference between 0.5 and 1 .25 of the mixing- 
length parameter (Castelli et al., 1997). Generally, it's very dif- 
ficult to state the location of the continuum across the wings of 
the Balmer lines, especially in Echelle spectra which are noto- 
riously severely affected by the blaze function. 

In computing the synthetic spectra we used the Kurucz's 
atomic and molecular line-lists (Kurucz 1992). The effects of 
several molecules were taken into account, including C2, CN, 
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Fig. 2. The effect of ff-enhancement, micro-turbulent velocity 
^, new vs old ODFs, and inclusion of 'predicted lines' are il- 
lustrated by sample spectra. To increase readability all spectra 
have been broadened to a rotational velocity of 20 km sec -1 . In 
each panel the bottom part highlight the difference between the 
plotted spectra (always in the sense solid - dotted lines). 
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Table 1. Comparison of the principal characteristics of publicly available libraries of synthetic spectra with resolution better 
than 5 A. N: number of combinations of atmospheric parameters and rotational velocities, new ODF: use of the new opacity 
distribution functions. TiO: inclusion of the TiO line-list. pred. lines: inclusion of the ' predicted lines'. The last three columns 
indicate if the spectra are provided for different values of enhancement of a-elements, micro-turbulence velocity and rotational 
velocity. 





atm. 


resol. 


N 




range 


of 




new 


TiO pred. 


variation 


of 




spectra 


power 




A(A) 


Teff(K) 


log 9 


[M/H] 


ODF 




lines 


[a/Fe] 




v rot 


Chavez et al. (1997) 


Kurucz 


250000 


711 


4 850 

5 400 


4 000 
8 000 


1.0 

5.0 


-1.0 
+0.5 


no 


no 


no 


no 


no 


no 


Zwitter et al. (2004) 


Kurucz 20000 61196 
11500, 8 500 


7653 
8 747 


3 500 
47500 


0.0 
5.0 


-3.0 
+0.5 


yes 


yes 


no 


yes 


yes 


yes 


Murphy & Meiksin (2004) 


Kurucz 


250000 


6410 


3000 
10 000 


5 250 
55 000 


0.0 
5.5 


-5.0 
+1.0 


no 


no 


yes 


no 


no 


no 


Martins et al. (2005) Kurucz 0.3A/pix 

Phoenix, Tlusty 


1654 


3000 
7000 


3 000 
55 000 


-0.5 
5.5 


-1.0 
+0.3 


no 


yes 


no 


no 


no 


no 


Rodriguez-Merino et al. (2005) 


Kurucz 


50 000 


1770 


850 
4 700 


3 000 
50 000 


0.0 
5.0 


-2.0 
+0.5 


no 


no 


no 


no 


no 


no 


Brott & Hauschildt (2005) 


Phoenix 


2 A/pix 


44000 


2 000 
10000 


2 700 
10 000 


-0.5 
5.5 


-3.5 
+0.5 




yes 




yes 


no 


no 


this library 


Kurucz 


20000 
11500 


51288 


2 500 
10 500 


3 500 
47500 


0.0 
5.0 


-2.5 
+0.5 


yes 


yes 


no 


yes 


yes 


yes 



8500, 2 000 
1, 10 A/pix 



Table 2. Range of atmospheric parameters explored by our synthetic spectral library. 



temperature (K) 


3500 < T eff < 47 500 


steps of 250 K for T sf! < 10 000 K; proportionally larger for higher T cft 


gravity (cgs units) 


0.0 < logg < 5.0 




metallicity 


-2.5 < [M/H] < +0.5 




rotation velocity (km sec -1 ) 


0< V rot <500 


0,2,5,10,15,20,30,40,50,75,100 for r cff < 6000 K 
0,10,20,30,40,50,75,100,150,200,250,300,400,500 for T eff > 6000 K 


enhancement 


[a/Fe] = 0.0, +0.4 


[a/Fe]=+0.4 for [M/H]=-0.5, -1.0, -1.5 


micro-turbulent velocity (km sec -1 ) 


£= 1,2,4 


1 and 4 km sec -1 for [a/Fe]=+0.4 


ODFs 


new 


old where new ones not yet available 


predicted lines 


excluded 





CO, CH, NH, SiH, SiO, MgH and OH. The TiO molecular lines 
(from Schwenke 1998) were included in the computation when 
r e ff < 5000 K. The synthetic spectra are calculated and made 
available both as fluxes for unit wavelength and as normalized 
intensities. The latter have a different meaning from the nor- 
malization to the continuum of the observational spectroscopy 
jargon, where the location of the continuum is taken to be rep- 
resented by the regions away from absorption lines or bands. 
In synthetic spectra, a continuum normalized spectrum is ob- 
tained by dividing the absolute fluxed spectrum by its calcu- 



lated continuum. Noticeable differences are essentially limited 
to the head of the Paschen and Balmer series of hydrogen, and 
to the coolest stars dominated by molecular bands. 



3. The library 

The spectra in our library covers from 2500 to 10500 A 
and were calculated at a resolving power R P =/l/A/l=500000 
and then degraded by Gaussian convolution to lower resolv- 
ing powers and properly re-sampled to Nyquist criterion (the 
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FWHM of the PSF being 2 pixels), to limit the data volume 
and therefore facilitate the distribution. 

The highest resolving power for our library is T? P =20000, 
which is typical of moderately high resolution Echelle obser- 
vations, as the spectral database of real stars of Munari and 
Tomasella (1999) and Marrese et al. (2003). The library is 
also provided at the resolving powers of the SLOAN survey 
R P =2 000 (York et al. 2000), of the RAVE survey fi P =8 500 
(Steinmetz 2003, Munari et al. 2005) and of the Gaia space 
mission tf P =ll 500 (Katz et al. 2004, Wilkinson et al. 2005). 
Furthermore, the present library is also made available at 
the uniform dispersions of 1 and 10 A/pix (again with the 
FWHM of the PSF being 2 pixels). The latter two are in- 
tended to be of assistance in classification of medium and low 
resolution grating spectra obtained in single dispersion mode 
(as with conventional B&C spectrographs or EFOSC-like im- 
ager/ spectrographs) . 

Each spectrum is provided for a range of rotational veloci- 
ties, as detailed in Table 1, spanning 14 values between and 
500 km sec -1 for stars hotter than 6000 K and 1 1 values be- 
tween and 100 km sec 1 for cooler ones. 

The explored combinations of parameters in Table 2 give 
rise to 5 1 228 spectra, each provided in six different combi- 
nations of resolving power and sampling, and two flux vari- 
eties (absolute units and normalized intensities). The total num- 
ber of individual spectra contained in this library is therefore 
51288x6x2=615456. The spectra are provided as gzipped 
ASCII files containing only the flux column, the common 
wavelength column being given separately only once for each 
of the six resolutions. The adopted scheme for the naming of 
files is outlined in Table 3. 

All spectra in the library can be directly ac- 
cessed and retrieved through the dedicated web page 
http://archives.pd.astro.it/2500-10500/ The version of the 
library at 1 A/pix is accessible also via ESA's web site 
http://gaia.esa.int/spectralib/, where browsing facilities based 
on Virtual Observatory tools are provided. A distribution via 
DVDs will be possible in special cases (to be arranged directly 
at munari@pd.astro.it). 

4. Test applications 

To evaluate the spectroscopic and photometric performances of 
our library we have carried out several tests. They support the 
fully satisfactory use of the library in spectroscopic applica- 
tions (our main goal), while caution applies on the photomet- 
ric use at the lower temperatures and/or shorter wavelengths. It 
is worth noticing that significant discrepancies affect the com- 
parison of available libraries with observational data over blue 
wavelengths and low temperatures. 

4.1. Eclipsing binaries 

A heavily demanding spectroscopic application of our library 
has been carried out by Siviero et al. (2004) and Marrese et 
al. (2005) in state-of-the-art investigation of double-lined so- 
lar type eclipsing binaries, as a source for (i) best matching 
templates in cross-correlation determination of radial velocities 



Table 3. Naming scheme for the spectra contained in the library 



example : T04000G45M05 V0 1 5K2 ANWN VR20F. asc 

characters meaning 

1-6 T04000 7; ff (K) 
7-9 G45 10 x log g (cgs units) 

10 M sign for metallicity : M= '-' , P= '+' 

11-12 05 10x[M/H] 
13-16 V01 5 rotational velocity (km sec -1 ) 
17-18 K2 microturbulent velocity (km sec -1 ) 

19 A [ff/Fe] enhancement: S=0.0, A=+0.4 

20-21 NW NW = new ODF models, 

OD = old ODF models 
22-23 NV no overshooting 

24-26 R20 R20 = resolving power 20 000; 

RVS = resolving power 11 500 (GAIA); 
RAV = resolving power 8500 (RAVE); 
SLN = resolving power 2000 (SLOAN); 
D01 = 1 A/pix dispersion; 
D10 = 10 A/pix dispersion 
27 F F=fluxed spectrum (erg cm -2 sec -1 A -1 ) 

N=normalized spectrum 



and (ii) reference grid in^ 2 analysis of atmospheric properties 
of binary components. They used the /?p=20000 version of the 
library. 

The availability of a wide and complete range of rotational 
velocities for all the spectra in our library always allowed to se- 
lect the best template even with binaries characterized by high 
rotations, a key ingredient to reach absolute precisions better 
than 1 km sec -1 in individual radial velocities even on modest 
S/N spectra. 

The derivation of atmospheric parameters via x 2 fitting to 
the library converged independently to the same values as de- 
rived geometrically by the orbital solution. The x 2 and orbital 
determinations provide the same gravities within 0.03 dex, the 
same rotational velocities within 4 km sec -1 and the same dif- 
ference in temperature between the two components within 
42 K (0.6%). The metallicity was indirectly derived by the or- 
bital solution via matching with theoretical isochrones tracks 
on the temperature-luminosity plane. The values derived in this 
way are confirmed within 0.06 dex by direct^ 2 analysis. Given 
the significant rotation exhibited by components of binaries 
(generally synchronous with orbital motion) the wide range of 
rotational velocities covered by our library proved once again 
to be a mandatory ingredient for accurate^ 2 results. 

The above results would have been significantly degraded 
if, instead of our one, it would have been used a library com- 
puted with inclusion of the PLs. To support this assertion we 
have re-run exactly the same x 2 fits of Siviero et al. (2004) 
and Marrese et al. (2005) using an identical grid of synthetic 
spectra, this time however computed on purpose with the in- 
clusion of PLs. In the case of V570 Per the use of spectra 
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V432 Aur HJD = 689.3270 * = 0.75 [M/H]= 0.6 




no PL 



i i i i i i i 

4700 4710 4720 X(A) 



V570 Per HJD-869.4234 $-0.62 [M/H]- + 0.02 




no PL - 



i i i i i i i i 

4700 4710 4720 X(A) 

Fig. 3. The observed spectra (solid thick lines) of two eclips- 
ing binaries are compared with synthetic ones computed with- 
out (solid thin line) and with (dotted line) inclusion of so 
called 'predicted lines'. The synthetic spectra are computed for 
the atmospheric parameters (T s g, logg) and co-rotational rota- 
tional velocity as derived from the accurate orbital solutions of 
Siviero et al. (2004, for V432 Aur) and Marrese et al. (2005, for 
V570 Per) and their chemical abundance analysis. The spectra 
of the two components are shifted at the appropriate radial ve- 
locity and combined together following the light-ratio at the 
these wavelengths as derived in the orbital solution. It is evi- 
dent from the comparison how the synthetic spectra computed 
without the 'predicted lines' offer a much better representation 
of the observed ones, and how the effect of the inclusion of 
'predicted lines' get worst with increasing metallicity. 



with PLs pushes the x 2 fit away from the geometrical solu- 
tion by Ar eff =109 K and A[M/H]=0.16, and by Ar eff =98 K, 
Alogg=0.10 and A[M/H]=0.11 for V432 Aur (average values 
for the components of the given binary). The differences are al- 
ways in the sense that inclusion of PLs pushes the x 2 fit away 
from the correct one by going toward cooler temperatures and 
lower metallicities. This is an expected, systematic trend. In 
fact, the inclusion of PLs causes the continuum in synthetic 
spectra to become fainter (cf. bottom panel of Figure 2), caus- 
ing a reduction in the contrast between strong absorption lines 
(those dominating the outcome of the x 2 fit) and the adjacent 
continuum. To re-establish the proper contrast, the x 2 fits a 
cooler temperature. Also fitting a lower metallicity helps in re- 
establish the proper contrast. In fact, PLs are mainly a forest of 



weak, optically thin lines that respond much faster to a small 
change in metallicity than the stronger absorption lines which 
are much closer to optically thick conditions than PLs. 

We observed also a reduction in the accuracy of the ra- 
dial velocities, which is however more difficult to quantify. In 
fact, it appears as a systematic shift (toward shorter or longer 
wavelengths) dependent on the actual wavelength interval and 
the orbital phase. This is caused by the beating in the cross- 
correlation between the pattern of PLs and the observed absorp- 
tion lines which is function of both A and radial velocity sep- 
aration between the binary components. The presence of these 
locally different systematic shifts become evident only when 
radial velocities from various wavelength intervals on the same 
spectrum are inter-compared (like those provided by individual 
orders in an Echelle spectrum). The overall effect is an increase 
of the scatter of individual measurements, different in amount 
from spectrum to spectrum. 

This concerns the analysis performed on wavelength re- 
gions dominated by weak PLs, as in the example on the bottom 
panel of Figure 2. There are however limited wavelength in- 
tervals where PLs much stronger than usual are stochastically 
clustered together. The result of ax 2 fit performed on any such 
interval would provide non-sense results. One such interval is 
illustrated in Figure 3, where the observed spectra of V432 Aur 
and V570 Per are taken from Siviero et al. (2004) and Marrese 
et al. (2005). While the fit with synthetic spectra computed 
without PLs provide - for the values of the atmospheric param- 
eters derived by the geometrical orbital solution - an almost 
perfect match to the observed spectra, the same spectra com- 
puted with inclusion of PLs bear no resemblance at all to the 
observed spectrum. It is evident that the photometric calibra- 
tion of a color-index that includes a narrow band centered on 
such a wavelength interval would provide non-sense results if 
carried out on synthetic libraries computed with PLs. 

4.2. RAVE spectra 

We are engaged in the data analysis of the RAVE spectro- 
scopic survey of the galactic Halo (Steinmetz 2003, Munari 
et al. 2005) with the 6dF fiber positioner and spectrograph at 
the UK-Schmidt at the Anglo- Australian Observatory. The sur- 
vey is carried out at high S/N over the 8420-8780 A interval at 
/?p=8 500 resolving power. 

Two low galactic latitude calibration fields rich in HD stars 
(1607m492 and 1716m425) have been observed a few times 
during the survey. Spectral classification of these stars are taken 
from the Michigan Project (Houk 1978). To test the perfor- 
mance of our library, we have^ 2 fitted the spectra of these HD 
stars against our library of synthetic spectra in its /? P =8 500 ver- 
sion. The results of the^- 2 fitting are presented in Figure 4. The 
right panels illustrate the accuracy of the achieved fits by over- 
plotting observed and synthetic spectra for the sample cases of 
a Sun analogue and of a A-type star. The T e g and log g values 
derived by the x 2 fit well match those expected from the G2V 
and A7IV classifications following the calibration by Straizys 
and Kuriliene (1981). For the whole sample of HD stars ob- 
served in the two calibration fields, the left panels of Figure 4 
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Fig. 4. The figure illustrates the use of our synthetic spectral library at resolving power /?p=8500 to classify RAVE spectra. 
The left panels show the results of x 2 fitting in gravity and temperature of spectra of HD stars from RAVE fields 1607m492 
and 1716m425, with spectral types from the Michigan Project (Houk 1978) and calibration into T e tf and logg by Straizys and 
Kuriliene (1981). The dotted lines represent margins for +3% in r e ff and ±0.25 dex in gravity. The error bars on each point 
represent the spread in the results ofx 2 fitting of three different spectra of the same object secured in different nights. On the right 
panels, examples of match between observed and synthetic spectra are provided for a couple of HD stars from the same RAVE 
fields. The results of the^ 2 fitting for the two spectra are given for T e g, logg, and metallicity. 



compare the T e g and log g values derived by the^- 2 fitting with 
those expected on the base of the spectral classification and 
the Straizys and Kuriliene (1981) calibration. The results in 
Figure 4 show an agreement within 3% for T e $ and better than 
0.25 dex in logg, the dispersion being almost entirely account- 
able for by the natural width in both T e s and log g of a spec- 
tral sub-type and a luminosity class in the MK classification 
scheme. 

4.3. The distance to the Pleiades 

The present library, in the 1 A/pix version, was used in dif- 
ferential broad-band photometric applications by Munari et al. 
(2004) while deriving a precise distance to the Pleiades via or- 
bital solution of the eclipsing binary member star HD 23642, 
which spectral type is A2V. 

Magnitudes in ten different photometric systems were 
found in literature for this star, and they were used to derive 
temperature, gravity and reddening of both components of the 
binary via calibration on our spectral libraries of the color in- 



dexes of each photometric system. The zero point of colors 
was fixed to 0.00 for a Vega model characterized by T e g =9400, 
log g=3.9, [M/H]=-0.5 and £=0.0 km sec 1 downloaded from 
Kurucz's web site. 

The photometric results on atmospheric and reddening pa- 
rameters so derived for both components of HD 23642 were 
confirmed by high resolution spectroscopy and the orbital so- 
lution. The successful photometric application of the library 
was made possible by the proximity in the parameter space of 
HD 23642 and the calibrator Vega. 

4.4. Broad-band photometric tests 

To test the broad-band photometric performances of our library 
we have investigated the reproduction of observed colors along 
the Main Sequence (MS). The results are shown in Figures 5 
and 6, where the r e tf and log g of stars along the main sequence 
were adopted from Straizys and Kuriliene (1981). We stopped 
at M2V as the MS faint end, because later types are charac- 
terized by r eff cooler than the 3500 K limit of our library (cf. 
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Fig. 5. Comparison of synthetic colors along the Main Fig. 6. The figure compares the behaviour of photometric col- 
Sequence truncated at M2V. The observed empirical Main ors alon 8 the Main Sequence for individual spectral types. The 
Sequence of Fitzgerald (1970) serves as a comparison. The same data of Fi 8 ure 5 are used - The ^opted band transmission 
adopted band transmission profiles are those of Buser (1978). P rofiles for *c and / c are those of Bessell (1990). 



Table 4). We included as a term of comparison all the other 
libraries from Table 1 for which data are actually accessible 
and are characterized by a wavelength coverage wide enough 
to cover at least three consecutive photometric bands. In this 
regard, a word of caution is necessary when considering the 
results from Martins et al. (2005) library, which red limit is 
placed at 7000 A where the V band has still 1 % transmission 
and does reach 0.0 only at 7400 A. The effect of missing this 
far red wing of the V transmission profile has effects on colors 
of cooler stars by making their B - V color marginally bluer. 

To assist in the interpretation of the result of the photomet- 
ric tests we have also recomputed the spectra of our library that 
define the MS, this time with inclusion of PLs (these spectra 
are available on the library web page). In Figure 5 and 6 by 
'NOVER 1221 pt' and 'ODFNEW 1221 pt' we mean the flux 
distributions that are predicted directly from the model atmo- 
sphere code, not to be confused with the spectra generated by 
the spectral synthesis code. 

Figure 5 illustrates the behaviour on classical 2-color dia- 
grams. Our library perform similarly well to the others on the 



U - B, B - V plane for colors bluer than 0.7 in B - V. For 
redder colors, the synthetic MS has U - B colors too blue by 
some hundredths of a magnitude. The inclusion of PLs in the 
computation brings the synthetic MS in close agreement with 
the observed one and very close to the ODFNEW 1221 pt flux 
distributions. The effect of adopting the ODFNEW models at- 
mospheres affects the lowest temperatures where the synthetic 
spectra (with and without PLs, as well as the 1221 pt flux dis- 
tributions) fails to match the very red colors observed in real 
stars, with the synthetic MS displaying a small upward hook. 

Figure 6 illustrates the data in a different format and ex- 
panded scale. The difference in the single color between real 
stars and synthetic spectra are plotted as function of the spec- 
tral type along the MS. This diagram highlights better the scat- 
ter between observed colors and those from synthetic spectra, 
particularly in U - B. It shows that, in all colors, significant de- 
partures of our library from observations are limited to spectral 
types cooler than K4V (r eff <4500 K). 

Having already proved the satisfactory spectroscopic per- 
formance of our library, we refrain from discussing any fur- 
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Table 4. Comparison between different calibrations for colors 
and temperatures of late type main sequence stars, a: Bessell 
(1990); b: Bessell (1995); c: Strayzis (1992); d: Bertone et al. 
(2004) using Kurucz spectra; <?: Bertone et al. (2004) using 
Phoenix spectra; /: Drilling and Landolt (2000); g: Straizys 
and Kuriliene (1981). 



Teff (K) 



Re — Ic color 
(«) (b) (c) (b) (d) (e) (/) (g) 



K5V 0.63 0.62 0.55 

MOV 0.85 0.91 0.89 

M2V 1.08 1.16 1.11 

M4V 1.36 1.63 1.42 



4350 4480 
3900 4010 
3600 
3130 



4630 4410 4395 
3990 3840 3900 
3720 3520 3550 
3450 3250 



ther its broad-band photometric performances and that of the 
other libraries listed in Table 1, which is clearly beyond the 
scope of this paper. Figures 5 and 6 should not be considered 
as the ultimate comparison of the photometric performances of 
the available synthetic libraries, even if we placed all libraries 
on the same reference system, using the same UBVRqIc band 
profiles and the same r e ff/log g conversion to spectral type for 
all of them. Even if different reconstruction of the UBVRqIc 
band profiles or r ef f/log g conversion would have been used, the 
spread among different synthetic libraries and with respect to 
observed colors would have remained. As a matter of the cau- 
tion that one has to bear in main in such comparisons, Table 4 
shows the widely different T e s calibration that have been pub- 
lished for the cooler spectral types along the MS considered in 
Figures 5 and 6. 
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Fig. 7. Figure la, to be published electronic only. 
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Fig. 9. Figure lc, to be published electronic only. 
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Abstract. To investigate the physical nature of the 'nucleated instability' of proto giant planets (Mizuno 1980 1, the stability of 
layers in static, radiative gas spheres is analysed on the basis of Baker's 1966 standard one-zone model. It is shown that stability 
. depends only upon the equations of state, the opacities and the local thermodynamic state in the layer. Stability and instability 

can therefore be expressed in the form of stability equations of state which are universal for a given composition. 

in . 

The stability equations of state are calculated for solar composition and are displayed in the domain -14 < lgp/[gcm 3 ] < 
0, 8.8 < lge/[ergg~'] < 17.7. These displays may be used to determine the one-zone stability of layers in stellar or planetary 
structure models by directly reading off the value of the stability equations for the thermodynamic state of these layers, spee- 
ch) i 1 ified by state quantities as density p, temperature T or specific internal energy e. Regions of instability in the (p, e)-plane are 
described and related to the underlying microphysical processes. Vibrational instability is found to be a common phenomenon 
at temperatures lower than the second He ionisation zone. The ^--mechanism is widespread under 'cool' conditions. 



CO 

Key words, giant planet formation - zr-mechanism - stability of gas spheres 



1. Introduction 

In the nucleated instability (also called core instability) hypothesis of giant planet formation, a critical mass for static core 
envelope protoplanets has been found. Mizuno ( 1980) determined the critical mass of the core to be about 12 M ffi (M e = 5.975 x 
10 27 g is the Earth mass), which is independent of the outer boundary conditions and therefore independent of the location in the 
solar nebula. This critical value for the core mass corresponds closely to the cores of today's giant planets. 

Although no hydrodynamical study has been available many workers conjectured that a collapse or rapid contraction will 
ensue after accumulating the critical mass. The main motivation for this article is to investigate the stability of the static envelope 
at the critical mass. With this aim the local, linear stability of static radiative gas spheres is investigated on the basis of Baker's 
( 1966 1 standard one-zone model. 



Phenomena similar to the ones described above for giant planet formation have been found in hydrodynamical models con- 
cerning star formation where protostellar cores explode (Tscharnuter 1987 Balluch [i988> . whereas earlier studies found quasi- 
steady collapse flows. The similarities in the (micro )physics, i.e., constitutive relations of protostellar cores and protogiant planets 
serve as a further motivation for this study. 



Fig. 1. Adiabatic exponent Y\. Y\ is plotted as a function of lg internal energy [ergg : ] and lg density [g cm 3 ]. 
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2. Baker's standard one-zone model 



In this section the one-zone model of Baker J1966I . originally used to study the Cepheid pulsation mechanism, will be briefly 
reviewed. The resulting stability criteria will be rewritten in terms of local state variables, local timescales and constitutive 
relations. 



Baker ( 1966 1 investigates the stability of thin layers in self-gravitating, spherical gas clouds with the following properties: 

- hydrostatic equilibrium, 

- thermal equilibrium, 

- energy transport by grey radiation diffusion. 

For the one-zone-model Baker obtains necessary conditions for dynamical, secular and vibrational (or pulsational) stability (Eqs. 
(34a, b, c) in Baker 1966). Using Baker's notation: 

M r mass internal to the radius r 

m mass of the zone 

ro unperturbed zone radius 

po unperturbed density in the zone 

7q unperturbed temperature in the zone 

L r o unperturbed luminosity 

£th thermal energy of the zone 

and with the definitions of the local cooling time (see Fig.^ 
_ E t h 

fco — T > 

and the local free-fall time 



(1) 



Tff 



3n 4?rrg 



Baker's K and ctq have the following form: 
it 1 
V8 T ff 
V32 1 T_ ff 

n 6 r co 

where « m(Po/po) has been used and 



0"o 



K = : 



e — mc 



(2) 

(3) 
(4) 

(5) 



is a thermodynamical quantity which is of order 1 and equal to 1 for nonreacting mixtures of classical perfect gases. The physical 
meaning of <x and K is clearly visible in the equations above. cr represents a frequency of the order one per free-fall time. K 
is proportional to the ratio of the free-fall time and the cooling time. Substituting into Baker's criteria, using thermodynamic 
identities and definitions of thermodynamic quantities, 



Ti = 



'ad 



( d\n_P \ 
[d\np) s 

I dhxT 
UlnP 



Xp 



Xt 



d\nP\ 
d\npj T ' 

d\nP 
d\nT 



dlnP 

d\\\K 
d\nT 



one obtains, after some pages of algebra, the conditions for stability given below: 

n 2 1 



-(3T] -4) > 



'ft' 



-r,v 



1 V a d 



1 - 3/4* p 



{k t - 4) + k p + 1 



' CO 1 ff 



4V ad - (Vad/C r + Kp) - — 



> 



> 



(6) 



(7) 



(8) 
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Table 1. Opacity sources. 
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Source 


77 [K] 


Yorke 1979, Yorke 1980a 


< 1700 a 


Kriigel 1971 


1700 < T < 5000 


Cox & Stewart 1969 


5000 < 



This is footnote a 



For a physical discussion of the stability criteria see Baker ( 1966) or Cox ( 1980 1. 

We observe that these criteria for dynamical, secular and vibrational stability, respectively, can be factorized into 

1 . a factor containing local timescales only, 

2. a factor containing only constitutive relations and their derivatives. 

The first factors, depending on only timescales, are positive by definition. The signs of the left hand sides of the inequalities 0, 

and © therefore depend exclusively on the second factors containing the constitutive relations. Since they depend only on 
state variables, the stability criteria themselves are functions of the thermodynamic state in the local zone. The one-zone stability 
can therefore be determined from a simple equation of state, given for example, as a function of density and temperature. 
Once the microphysics, i.e. the thermodynamics and opacities (see Table[Q, are specified (in practice by specifying a chemical 
composition) the one-zone stability can be inferred if the thermodynamic state is specified. The zone - or in other words the layer 
- will be stable or unstable in whatever object it is imbedded as long as it satisfies the one-zone-model assumptions. Only the 
specific growth rates (depending upon the time scales) will be different for layers in different objects. 

We will now write down the sign (and therefore stability) determining parts of the left-hand sides of the inequalities Q 
and and thereby obtain stability equations of state. 

The sign determining part of inequality is 3r"i - 4 and it reduces to the criterion for dynamical stability 

rt>l. (9) 

Stability of the thermodynamical equilibrium demands 

* p >0, c v >0, (10) 
and 

Xt >0 (11) 
holds for a wide range of physical situations. With 

r 3 - 1 = 4— > ( 12 ) 

pT c v 

Ft =* p +*r(r 3 -l) > (13) 

V ad = T - 1 ^- > (14) 
1 1 

we find the sign determining terms in inequalities @ and (jSJi respectively and obtain the following form of the criteria for 
dynamical, secular and vibrational stability, respectively: 

3r 1 -4=:5 dyn > (15) 

1 - 3/4^ p 

-{k t -A) + k p + 1 =: S sec > (16) 

Xt 

4 

4V ad - (V ad * r + k p ) - — =: S vib > . (17) 



3r 



The constitutive relations are to be evaluated for the unperturbed thermodynamic state (say (po, To)) of the zone. We see that 
the one-zone stability of the layer depends only on the constitutive relations Fi, V ad , x T , x p , k p , k t . These depend only on the 
unperturbed thermodynamical state of the layer. Therefore the above relations define the one-zone-stability equations of state 
S dyn , Ssec and Svib- See Fig.|2]for a picture of 5 v ib- Regions of secular instability are listed in Table 1. 
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Fig. 2. Vibrational stability equation of state 5 , v ib(lg e, lgp). > means vibrational stability. 

3. Conclusions 

1. The conditions for the stability of static, radiative layers in gas spheres, as described by Baker's ( 1966 1 standard one-zone 
model, can be expressed as stability equations of state. These stability equations of state depend only on the local thermody- 
namic state of the layer. 

2. If the constitutive relations - equations of state and Rosseland mean opacities - are specified, the stability equations of state 
can be evaluated without specifying properties of the layer. 

3. For solar composition gas the A:-mechanism is working in the regions of the ice and dust features in the opacities, the H2 
dissociation and the combined H, first He ionization zone, as indicated by vibrational instability. These regions of instability 
are much larger in extent and degree of instability than the second He ionization zone that drives the Cepheid pulsations. 
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